function [criterion diff weight] = evaluation(empmom,simmom,covar,momlist)

MOM = zeros(2,length(momlist));
weight = zeros(length(momlist));
for j = 1 : length(momlist)
 MOM(1,j) = empmom(momlist(j));
 MOM(2,j) = simmom(momlist(j));
 for jj = 1 : length(momlist)
  weight(j,jj) = covar(momlist(j),momlist(jj));
 end
end

diff      = MOM(1,:) - MOM(2,:);
criterion = diff*weight^(-1)*diff';